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ABSTRACT 

We reanalyze r-process nucleosynthesis in the neutron-rich ejecta from a prompt 
supernova explosion of a low-mass (11 M©) progenitor. Although it has not yet been 
established that a prompt explosion can occur, it is not yet ruled out as a possibility 
for low-mass supernova progenitors. Moreover, there is mounting evidence that a new 
r-process site may be required. Hence, we assume that a prompt explosion can occur 
and make a study of r-process nucleosynthesis in the supernova ejecta. To achieve a 
prompt explosion we have performed a general relativistic hydrodynamic simulation of 
adiabatic collapse and bounce using a relativistic nuclear-matter equation of state. The 
electron fraction 1^ during the collapse was fixed at the initial-model value. The size of 
the inner collapsing core was then large enough to enable a prompt explosion to occur 
in the hydrodynamical calculation. Adopting the calculated trajectories of promptly 
ejected material, we explicitly computed the burst of neutronization due to electron 
captures on free protons in the photodissociated ejecta after the passage of the shock. 
The thermal and compositional evolution of the resulting neutron-rich ejecta originating 
from near the surface of proto-neutron star was obtained. These were used in nuclear 
reaction network calculations to evaluate the products of r-process nucleosynthesis. We 
find that, unlike earlier studies of nucleosynthesis in prompt supernovae, the amount 
of r-process material ejected per supernova is quite consistent with observed galactic 
r-process abundances. Furthermore, the computed r-process abundances are in good 
agreement with Solar abundances of r-process elements for A> 100. This suggests that 
prompt supernovae are still a viable r-process site. Such events may be responsible for 
the abundances of the heaviest r-process nuclei. 
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1. Introduction 

The astrophysical site for the origin of elements 
by rapid neutron capture (the r-process) has been a 
mystery for some time. Numerous sites have been 
proposed as possible candidate environments (Hille- 
brandt 1978; Mathews & Ward 1985; Mathews & 
Cowan 1990; Meyer 1994). At the present time 
thei most popular among proposed sites involves the 
neiitrino-hcated ejecta from a nascent neutron star 



(Woosley & Hoffman 1992; 


Meyer et al. 1992; 


Woosley 


3t al. 1994 


; Otsuki et al. 2000; 


Sumiyoshi et al. 200C|; 


lerasawa et al. 2001 


). In Woosley et al. (1994) it was 



demonstrated that the Solar r-process abundances 
were well reproduced in material which has been ab- 
lated from the proto-neutron star in neutrino-driven 
winds. 

There are, however, a few serious problems with 
this r-process paradigm. For one, elements with 
A ~ 90 are significantly overproduced by over a factor 
of ~ 100. Second, the requisite high entropy {S/k > 
400) in the supernova simulations has not been du- 
plicated by other independent theoretical studies (cf. 



Witti, Janka, & Takahashi 1994; Takahashi, Witti, & 



Janka 1994; Qian & Woosley 1990). A third prob- 



lem (Meyer, McLaughlin, & Fuller 1998) is the possi- 



bility that neutrino-nucleus interactions could drasti- 
cally alter the neutron to seed ratio rendering a sat- 
isfactory r-process more difficult. 

On the other hand, other viable sites have been 
demonstrated to also account for the Solar r-process 
abundance pattern. For example, it has been s hown 
( Freiburghaus, Rosswog, fc Thielemann 1999| ) that 



neutron-star mergers can produce the abundances of 
nuclei with A > 130. Collapsing 0-Ne-Mg cores re- 
sulting from progenitor stars of low mass ~ IOMq 
also remain ( Wheeler, Cowan, fc Hillebrandt 1998 ) as 



a promising site for the r-process. The present work, 
however, is primarily concerned with collapsing iron 
cores of ^ llAf© progenitor stars. 

In several early papers (e.g. Hillebrandt et al. 1976; 
Hillebrandt 1982; Hillebrandt et al. 1984) studies were 
made of r-process nucleosynthesis in the material hy- 
drodynamically ejected during a prompt supernova 
explosion. In Hillebrandt et al. (1976) the r-process 
was thought to occur in the inner ~ 0.43 Mq of 
ejected material in which the neutron to proton ratio 
was taken to range from n/p ~ 1 to 7. The prob- 
lem with these earlier studies, however, is that it ap- 
peared that too much material was ejected to make a 



reasonable accounting for the galactic abundances of 
r-process nuclides. To account for the current galac- 
tic abundance of r-process material requires that only 
about 10""* Mq of r-process elements be ejected per 
supernova (cf. Mathews & Cowan 1990; Woosley et 
al. 1994). This is much less than the (~ 0.4 Mq) 
ejected in the Hillebrandt et al. (1976) model. Fur- 
thermore, it gradually became clear (p3aron & Coop-I 



erstcin 1990, Bcthe 1990, Suzuki 1994) that models 



of prompt supernovae could only be made to explode 
for a limited range of small core-masses, high lepton 
fractions, and equation of state parameters, and may 
not occur at all. Because of these dilemmas, prompt 
supernovae have fallen out of favor as an r-process 
paradigm in recent years. 



However, recent studies ( Hcger et al. 2001 ) of pro- 
genitor cores b ased upon improved slower electron 
capture rates ( Langanke fc Martinez-Pincdo 2000 ) 
have shown that smaller, cooler, iron cores with larger 
lepton fractions are a likely outcome of supernova 
progenitor-star evolution. Although it has not yet 
been established that a prompt explosion can occur, 
the newer electron capture rates may suggest that 
they are possible. 

The purpose of the present work is, therefore, to 
revisit the question of whether prompt supernovae 
necessarily overproduce the abundances of r-process 
elements. We assume that a prompt explosion can 
occur under the plausible conditions described below. 
We then examine the associated r-process nucleosyn- 
thesis. Our purpose is not to provide a fully detailed 
prompt explosion model, which will take some time to 
complete. Instead, our purpose here is merely to de- 
velop a basic prompt explosion model with enough 
content to examine the question of the implied r- 
process yield and relative abundances. Therefore, we 
begin with an adiabatic collapse simulation with the 
electron fractions (Fg) of the various zones fixed at 
the initial-model values. This approximates the effect 
of the new diminished electron capture rates. This 
model produces a prompt explosion, both because the 
core bounce is energetic, and because the higher 
leads to a small outer iron core which does not have 
enough material to dissipate the shock by photodis- 
sociation. It is adequate for our purpose. 

A key part of the r-process calculation of the 
present work is that we perform post processing of 
the compositional change of the ejecta relevant to de- 
termining the amount of r-process material ejected. 
First, we consider the burst of neutronization due to 
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electron capture on free protons after the shock heat- 
ing and photodissociation of the ejecta. We show that 
the neutronization is rapid. Moreover, the amount 
of neutronized ejecta is much less than was assumed 
in earlier work. Consequently, the correct amount of 
neutron rich material is ejected to account for galactic 



abi ndances. Our network calculations of the ensuing 



r-process nucleosynthesis makes a reasonable repro- 
duction of the Solar r-process abundance distribution 



forlnucici witti A ^ iUU. i'tus suggests that prompt 



supernovae remain as a viable candidate site tor tUe 
r-process synthesis of heavy elements. 

2. Prompt Supernova Mechanism 

For some time it has been debated in the literature 
as to whether collapsing iron cores of supernova pro- 
genitors can explode via a prompt mechanism ( Bethe 
1990, Suzuki 1994). It has not yet been established 



that a prompt explosion can occur. The question has 
been whether the core bounce itself can be sufficiently 
energetic to eject the outer layers of the star. It is 
always difficult to develop a supernova model which 



actually explodes (Bcthc 199C, Suzuki 1994) and it 
has been clear for some time that if a prompt explo- 
sion is to occur it must involve a low-mass 10 — 12 
Mq) progenitor star and a low-m ass collapsing core 
(e.g. Baron fc Cooperstein 1990 ). Whether or not 



an explosion can occur depends upon a difference be- 
tween two large numbers: the energy in the outward 
moving shock and the energy lost by the shock due to 
the photodissociation of nuclei. At present there are 
too many uncertainties in the input supernova physics 
to definitively exclude a prompt explosion mechanism 
for at least a limited range of core masses and equa- 
tion of state parameters. Nevertheless, it is clear what 
ingredients are necessary. Two key factors which we 
invoke here to obtain a prompt explosion are a small 
total iron core mass and a large Ye as a result of slower 
electron capture rates. The important roles of these 
we now outline in a brief summary of the explosion 
mechanism. 



\_s the electron-degenerate iron core reaches its 



Chandrasekhar mass, the collapse divides into an in- 
ner homologously collapsing core and an outer core 
collapsing more slowly. A prompt explosion requires 
a small outer iron core so that the outward moving 
shock from the bounce of the inner core is not dissi- 
pated by photodissociatiori of the outer core. A small 



outlnr core can be a,chieved by minimizing the total 



Chandrasekhar mass of the iron core, while maxi- 
mizing the mass of the inner homologous core. The 
Chandrasekhar mass is smallest for the cooler cores 
of the lowest-mass supernova progenitors. A suffi- 
ciently small, cool, iron core {Mcore ~ 1-3 Mq) results 
from the 11 Mq progenitor star model of (Woosley & 



Weaver 1995) employed in the present work. 



A large inner homologous core can result if the lep- 
ton fraction is large (Takahara & Sato 1984, Yamada 



& Sato 1994). This can happen if the electron capture 



rates are slower than the usually employed rates of 
Fuller, Fowler, & Newman (1980;1982ab;1985) [here- 
after FFN]. Electron capture rates are important be- 
cause degenerate electrons account for a large fraction 
of the pressure support of the inner core. A larger 
therefore implies more pressure support, and a larger 
inner core. Ultimately, during the collapse, the neu- 
trinos become trapped. The electron capture rates 
are then suppressed by the high neutrino Fermi en- 
ergy which limits the available phase space. After- 
ward, the lepton fraction remains fixed. 

The collapse of the inner core slows as the central 
density exceeds nuclear density and an outward mov- 
ing shock is produced. A somewhat soft equation of 
state helps the explosion ( Baron fc Cooperstein 1990 ), 
but it is not as important as keeping the lepton frac- 
tion large. 

If the outer part of the core is small enough, it 
will not completely absorb the outward going shock 
by the photodissociation of iron-group nuclei. The 
shock is then not dissipated. It eventually reaches 
the surface of the outer iron core and an explosion 
can ensue. Therefore, this combination of conditions 
(i.e. a large lepton fraction and a small total iron core 
mass) enables energetic, prompt explosions and the 
ejection of neutronized material from deep inside the 
star. 

If the standard FFN electron capture rates are 
employed even smaller iron cores (M w 1.2 Mq) 
cannot explode promptly ( Rampp fc Janka 200d|, 



Liebendorfer et al. 2001; cf. Baron fc Cooperstein] 



199C). With these rates, the lepton fraction is too 
small and so is the inner core. However, with slower 
electron capture rates, a prompt explosion is possi- 
ble. The determination of electron capture rates in 
supernova cores has, however, suffered from large un- 
certainties in the associated nuclear physics. More- 
over, recent studies (e.g. Langanke fc Martinez 



Pinedo 2000 ) of electron capture on neutron- rich nu- 



clei have reported rates which are on average an order 
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of magnitude smaller than the previous standard rates 
of FFN. Indeed, a recent study ( |Heger et al. 200l| ) 
of the effect of these rates on 15-40 Mq progenitor 
stars found significantly larger and cooler, smaller, 
cores. Electron capture by free protons in a collapsing 
core can also be suppressed if the free proton fraction 
is small. A small proton fraction can result from a 
large symmetry energy in the equation of state. This 



i s t|ie case for the one wc adopt in the current study 



dShcn ct al. 1998a| , |Shcn ct al. 1998bD . 

Clearly, further numerical simulations of the for- 
mation and collapse of the small iron cores associ- 
ated with improved electron capture rates and a good 
equation of state are called for. This, however, is be- 
yo nd the scope of the present work. Nevertheless. 



thc l model described here is a plausible version of the 



required physics in order to study the r-process ele- 
ments produced. 

3. Explosion Model 

In the current study, we focus on the consequences 
for r-process nucleosynthesis of a prompt supernova 
explosion. In order to investigate mass ejection in 
a prompt explosion an adiabatic hydrodynamic col- 
lapse calculation was performed. As an extreme case, 
we assume that the electron captures are sufficiently 
slower than the standard FFN rates that the elec- 
tron fraction remains fixed during the collapse. This 
maximizes the shock energy, and therefore, the ex- 
plosion energy, and brings the iron core to a prompt 
explosion. We use this hydrodynamical calculation to 
evaluate the amount of neutron-rich material ejected 
and the thermodynamic conditions during nucleosyn- 
thesis. 

Our ultimate goal is a full calculation of the neutrino- 
radiation hydrodynamics. It is, however, computa- 
tionally difficult to follow this evolution through the 
long time scale over which the nucleosynthesis occurs. 
The work described below is a necessary first prelim- 
inary step toward that goal. 

3.1. Equation of State 

One important ingredient for a calculation of a 
prompt explosion is the nuclear equation of state. 
Fortunately, recent information on the physics of un- 
stable nuclei is helping to clarify both supernova ex- 
plosions and the r-process. One can, for example, 
now use radioactive nuclear beam facilities to probe 
both the neutron-rich matter of supcrnovac and the 



nuclear data needed for r-process nucleosynthesis cal- 
culations. Having this experimental information on 
unstable nuclei with large asymmetry, Sumiyoshi et 
al. (1993; 1995a; 1995b) and Hirata et al. (1997) 
have extensively studied nuclear structure and nu- 
clear matter within a relativistic many-body frame- 
work. Based upon this work a relativistic EOS ta- 



ble for application to supernova simulations ( 3hen et 
al. 1998a, phcn ct al. 1998b ) has recently been devel- 



oped. This relativistic EOS table enables one to per- 
form full simulations of supernovae from the initial 
gravitational core collapse to the ejection and cooling 
of material from the newly formed neutron star. It 
has successfully been applied to studies of neutrino- 
driven winds from proto- neutron stars (Sumiyoshi et 



al. 2000) and numerica l simulations of core collapse 
( Sumiyoshi et al. 2001 ). Although this equation of 
state is not particularly soft (incompressibility =281 
MeV), it does, have a large symmetry energy (36.9 
MeV). Since nuclei with large asymmetries were em- 
ployed in the development of this equation of state, 
one expects that the symmetry energy is better de- 
termined. The large deduced value reduces the free 
proton fraction. This can raise the for the core 
since there are fewer protons available for electron 
capture. 

3.2. Hydrodynamics 

The prompt supernova calculations described herein 
are based upon the general-relativistic spherically- 
symmetric hydrodynamics code of Yamada (1997). It 
makes use of a fully implicit Lagrangian numerical 
scheme. It has been designed for the study of super- 
nova explosions and treats both the hydrodynamics 
and the neutrino transport. 

Although a Boltzmann solver for the neutrino 
transport has recently been implemented (Yamada, 
Janka, & Suzuki 1999) and detailed simulations of 
gravitational core collapse are now being undertaken, 
in the current study we have performed calculations of 
core collapse using only pure hydrodynamics without 
neutrino transport. Since we are primarily interested 
here in the dynamics of collapse and prompt ejection, 
we can as a first approximation ignore the late time 
neutrino diffusion and heating which is not manifest 
till later times in the explosion. Indeed, the present 
work is mostly concerned with matter ejected within 
^ 0.2 sec after core bounce and largely before the ar- 
rival of the delayed neutrinos. Obviously, this is an 
approximation which we will address in a subsequent 
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work. 

The present study, however, should be adequate 
to analyze the physics of the hydrodynamics of core 
collapse, bounce and a possible explosion with our 
adopted EOS. In particular, our primary purpose here 
is to evaluate the amount of ejected material which 
contributes to the r-process. For this particular pur- 
pose, our neglect of neutrino radiative transport is 
probably adequate. 

For the same reasons, we assume that the collapse 
is adiabatic collapse. That is wc neglect the heat- 
ing and cooling via neutrinos. Thus, the entropy per 
baryon of each mass element remains constant dur- 
ing the collapse except for during the passage of the 
shock. 

Altogether, these approximations tend to maxi- 
mize the energy of the prompt explosion and provide 
an upper bound on the amount of ejecta. Thus, our 
estimates of ejected r-process material can perhaps 
be taken as an upper limit which is what we seek in 
the present study. 

3.3. Presupernova Model 

We adopt the presupernova model of a 11 Mq pro- 
genitor as computed by Woosley & Weaver (1995). 
This is the best case for a prompt explosion as it 
provides the smallest iron core {Mcore = 1-32 Mq). 
It is also the smallest progenitor mass among their 
models ranging from IIM© to AOMq. Both of these 
conditions arc favorable for a prompt explosion. For 
our hydrodynamic simulations we are only concerned 
with the central 1.62 Mq, which contains the iron 
core and outer layers. Wc adopt densities, electron 
fractions, temperatures and radii as a function of the 
baryon-mass mesh from the progenitor model. Other 
quantities are derived using our adopted relativistic 
EOS table. 

The simulations utilize 200 baryon-mass mesh points. 
This provides enough resolution to analyze the r- 
process in the ejecta. The inner 50 grid points span 
a range of baryon mass from up to 1.0 Mq. The 
outer 150 grid points span the remaining 0.62 Mq. 
The grid size is taken to be uniform in each of these 
two regions. We have checked our results in simula- 
tions using different mass meshes and find that they 
do not appreciably alter the nucleosynthesis results. 

For the outer boundary conditions, we have fixed 
quantities at the values obtained at the cut position 
in the 11 Mq progenitor model. 



4. Results 

4.1. Collapse and Shock Ejection 

The central core of the initial presupernova model 
is already marginally gravitationally unstable. There- 
fore, when we start the numerical calculation with our 
adopted initial conditions, the inner core starts col- 
lapsing immediately. Because we fix Yg as described 
in the next subsection, the bounced inner core size is 
large enough for the shock to realize a prompt explo- 
sion with matter ejection. 

Figure 1 shows the evolution of various mass shells 
as a function of time during the collapse. Here we 
see that a well formed outward going shock devel- 
ops within ~ 100 ms after core bounce. This short 
time scale justifies our neglect of neutrino transport 
in this calculation, since delayed neutrinos would not 
be expected to affect the explosion till several hun- 
dred milliseconds after core bounce. 

We have extracted trajectories from the numerical 
prompt explosion results for use in the subsequent cal- 
culations of electron capture on protons and r-process 
nucleosynthesis. The trajectories with grid numbers 
ranging from 90 to 111 out of 200 grid points are re- 
sponsible for the r-process. These are used for the nu- 
cleosynthesis calculations. The main featiires of these 
trajectories are summarized in Figures 2a-2d. 

From Figures 2a and 2b we see that these mass 
shells are compressed during the bounce to densities 
in excess of 10^^ g cm"'^ and are heated to tem- 
peratures well above ~ 1 MeV due to the passage 
of the shock. This heating is more than sufficient 
to completely photodissociate this material into free 
neutrons and protons. Electron capture on the free 
protons then leads to a burst of neutronization. 

Figure 2c shows the entropy per baryon in the out- 
going ejecta. The obtained entropy {S/k ~ 10) is 
much less than that obtained in supernova wind mod- 
els (e.g. Woosley et al. 1994), but it is still sufficiently 
high that heating by the subsequent r-process nucle- 
osynthesis (Sato 1974; Hillebrandt et al. 1976) does 
not appreciably affect the hydrodynamic evolution of 
the ejecta. To see this, consider that the total Q-value 
of a typical beta-decay chain back to stability (e.g. 
i30cd isorjg) jg j^bout 20 MeV (i.e. - 3 x lO''^ erg 
per r-process isotope produced). The ejected mass is 
0.025 Mq or ~ 3 X 10^^ nuclei (assuming A ~ 100). 
Therefore, the total energy released is only 9 x 10^* 
erg, which is only ~ 1% of the explosion energy. In 
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earlier studies of r-process heating (Sato 1974; Hille- 
brandt et al. 1976) the entropy was much smaller, 
0.1 — 1, so that heating by the r-process was much 
more important for the dynamics of the ejecta. 

4.2. Electron Fraction in the Ejecta 

During the core collapse, we assume complete sup- 
pression of electron capture on neutron-rich nuclei 
and free protons in the supernova core. We fix 
to the original values of the progenitor core during 
the hydrodynamical calculation of the core collapse. 
(We calculate, however, compositional change after 
the shock passage separately, to be discussed below.) 
This assumption overestimates the electron fraction 
(or lepton fraction after the neutrino trapping) dur- 
ing the collapse. It provides us with a large inner 
core, and hence, enough shock energy for a successful 
prompt explosion. We expect that this approxima- 
tion might lead to an overestimate of the amount of 
ejected r-process material. If that is the case, then it 
is adequate for our purpose, which is to demonstrate 
that not too much r-process material is ejected. 

After the core bounce, the shock passes through 
the outer part of the iron core and dissociates iron- 
group nuclei into free neutrons and protons. Near 
to the remnant, the density reaches ~ 10^^ g cm~'^. 
Here, electron capture on free protons occurs abruptly 
because of the high ambient electron degeneracy. The 
material becomes immediately transparent to neutri- 
nos due to the dissociation of the nuclei which had 
been the cause of neutrino trapping. Neutrino ab- 
sorption, which is the reverse of the electron capture, 
is negligible. This process does not depend upon the 
details of neutrino transport. It corresponds to the 
well known neutronization burst in the flux of su- 
pernova neutrinos. It causes the material to quickly 
neutronize right after the shock passes. The result- 
ing neutron-rich region can be generally seen in many 
supernova simulations (e.g. Suzuki 1994). The in- 



ner regions near the nascent neutron star tend to be 
neutron-rich. It is this electron capture process in the 
ejecta that produces the neutron excess necessary to 
achieve r-process conditions. 

To treat this neutronization, we have calculated 
the time evolution of the electron fraction due to 
electron capture on free protons by postprocessing the 
results of the hydrodynamical calculation. This elec- 
tron capture drastically reduces the electron fraction 
and then ceases quickly due to the rapid expansion of 
the ejecta. The degree of neutronization depends on 



the compression of matter at the bounce. The gen- 
eral tendency is that material deeper in the ejecta has 
a higher electron degeneracy and a smaller electron 
fraction. The resulting electron fractions in various 
ejected mass shells then become initial conditions for 
our subsequent r-process nucleosynthesis calculations 
(see Figure 2d). 

We should point out that this explicit determina- 
tion of the neutronization of the ejecta is an improve- 
ment over the previous hydrodynamic calculations of 
Hillebrandt et al. (1976). In that work a continuous 
distribution in Yg between Yg = 0.5 at the base of the 
silicon-burning shell and Y^ = 0.125 at the remnant 
neutron star was simply assumed for the r-process 
initial conditions. As we shall see, our calculations 
indicate that a much smaller region of ncutronized 
material actually emerges from a prompt explosion. 

Figure 3 summarizes the electron fraction of mate- 
rial at the beginning of nucleosynthesis as a function 
of baryon mass coordinate in the ejecta. The dotted 
line displays the initial value of Kg before the neu- 
tronization. As we shall see the r-process occurs only 
in regions with Yg ~ 0.15 to 0.25. In this calculation 
that corresponds to only 0.025 Mq of material. This is 
much less than in the earlier estimates of Hillebrandt 
et al (1976). The difference between the present re- 
sults and their calculation can be traced to their sim- 
plified treatment of the neutronization of the ejecta. 
We stress again that we have performed hydrodynam- 
ical calculations to find the ejected part of material 
and have explicitly calculated the electron fractions of 
each mass elements due to electron captures on free 
protons in order to determine both the amount and 
neutron-richness of ejecta. 

The grid point 90 {Mb = 1.163M0 at zone center 
in baryon mass coordinate) is the inner-most mass 
shell of the ejected material. This trajectory has the 
smallest electron fraction {Y^ — 0.16 in Figure 3). It 
contributes the most to the r-process elements around 
the A=195 abundance peak. The grid points from 90 
to 95 {Mb = 1.184M0), leading to Yg <, 0.25 con- 
tribute to the r-process nucleosynthesis from A=130 
to A=190. The grid point 111 {Mb = 1.250Mq) has 
an electron fraction 0.45. This trajectory only pro- 
duces elements up to A = 70 and does not contribute 
to the r-process. Accordingly, we do not add outer 
mass shells (grid points > 112) in the following nucle- 
osynthesis calculations as they do not contribute to 
the abundances of r-process elements. 
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5. r-Process Nucleosynthesis 



As the neurronized material is ejecxed it quickly 



experiences charged-particle and neutron-caprure re- 
actions to produce seed material followed by neu- 
tron captures and beta decays along the r-process 
path. The r-process nucleosynthesis calculations in 
the present work are based upon the nuclear reaction 
network as described in Terasawa et al. (2001). In 
that work it was shown that both neutron-rich light- 
mass nuclei as well as heavy nuclei can play an im- 
portant role in the production of both seed nuclei and 
r-process elements. The path through light nuclei is 
most important in very neutron rich conditions such 
as in the high entropy neutrino-driven wind or in the 
most neutron-rich ejecta studied here. We thus make 
use of an extended nuclear reaction network which in- 
cludes all relevant nuclear reactions from protons to 
heavy actinides (Terasawa et al. 2001). The nuclear 



reaction network includes nuclei from beta-stability 
to the neutron-drip line. For 10 < Z < 94, we have 
used the network of Meyer et al. (1992), which in- 
cludes about 3000 nuclear species. We have extended 
this network to include lighter nuclei as well as al- 



most all charged-particle reactions for A< 28 (Kajino 
& I poyd 1990| ; |Orito et al. 1997| ) plus all {a, n) reac- 
tions up to Z— 36. For the nucleosynthesis calculation 
we have checked the results with and without cap- 



ture on nuclei ( Meyer, McLaughlin, fc Fuller 1 



The results do not depend upon this effect, presum- 
ably because material is so quickly removed from the 
neutron star surface. 

5.1. r-Process Yields 

The ejected mass shells of Figure 1 were evolved 
with our r-process network. Figure 4 shows the 
results of several trajectories labeled by their La- 
grangian zone numbers. We see that material ejected 
at the bottom contribute to the heaviest A = 195 
peak, while those shells which had higher Ye only con- 
tribute to the A = 130 peak and below. Shells with 
Ye > 0.25 only produce iron group and elements up 
to A = 100. They do not contribute to the synthesis 
of r-process elements. 

The nucleosynthesis yields of these different mass 
shells were summed to produce the final r-process 
abundance curve shown in Figure 5. Here, we see an 
adequate reproduction of the Solar r-process abun- 
dances. The dips above and below the A « 130 
peak in our computed abundance curve may be an 



indication (cf. [Meyer et al. 1992| ) that the N = 82 
shell closure is too strong in our mass formula (Hilf et 
al. 1976). Other than that, the abundances of nuclei 
heavier than A 100 are nicely reproduced. 

The fact that elements with A < 100 are under- 
produced suggests an interesting possibility. These 
are precisely the nuclei overproduced in the neutrino- 
driven wind scenario. They are also the elements 
which seem to deviate from the Solar r-process abun- 
dances in metal-poor halo stars (Ryan, Norris, & 
Beers 1996; Sneden et al. 1996; 1998; 2000). The 
suggestion is therefore that two different environ- 
ments may contribute to the r-process: one being 
the prompt explosion of low-mass SNe progenitors 
under study here. Such environments may produce 
the heaviest r-process nuclei; the other environment 
could be the neutrino-driven winds in more massive 
progenitors. These events might produce the lighter 
A < 100 nuclides. 

This is consistent with recent meteoritic evidence 
(Qian 2000; Qian & Wasserburg 2000) that two dif- 
ferent r-process environments preceded the formation 
of the Solar System. It is also consistent with galac- 
tic evolution evidence. Ishimaru & Wanajo (1999) 
have demonstrated that the observed large dispersion 
of r-process elements in metal-poor halo stars can be 
explained if the r-process occurred in either low-mass 
(8-10 Mq) or high mass (~ 30 Mq) supernova pro- 
genitors. Here we favor low-mass progenitors, but 
perhaps both environments contribute. 

5.2. Galactic Contribution 



Regarding the r-process yields implied by the present 
study, we see from Figure 3 that w 0.025 M© of 
r-process material is ejected in our 11 Mq model. 
This corresponds to a mass fraction in the ejecta of 
Xr{m = 11) « 3 X 10-3. However, since the prompt 
mechanism can only work in low mass progenitors, 
only supernovae in a narrow mass range, say Am ^ 1 
Mq probably contribute to the galactic abundances. 
Hence, to estimate the average yield per supernova, 
we must correct for the fraction of supernovae which 
can contribute. Assuming a standard Salpeter initial 
mass function (l){m) oc to^^ '^, then the total average 
mass fraction of ejected r-process material per super- 
nova {Xr)sN becomes: 



r SN 



Xr{m — 11) X (to — mrem)</'(ll)ATO 

J^{m - mrem)<l>{m)dm 



8 



0.3Js:^(m = ll)Am 

11 Mq 



7.5 X 10" 



,(1) 



where in the above we have ignored the mass of the 
remnant compared to the mass of ejecta. Now assum- 
ing that the total SNII rate has been constant at 10~^ 
yr~^ over the past 10^° yr, then we would estimate a 
current mass in r-process elements of the Galaxy of 
j^Gai 7 5 X -j^^^ Assuming a total baryonic 
mass for the Galaxy of 10^^ Mq, then the present 
mass fraction in r elements would be ^ 7.5 x 10"'' in 
reasonable agreement with the observed Solar mass 
fraction of ~ 10"'^. 

6. Conclusion 

We conclude that the prompt mechanism in low- 
mass supernova progenitors remains as a viable model 
for the r-process. Clearly, one must still establish 
whether prompt explosions of the sort described here 
actually occur when neutrino transport is included. 
Nevertheless, our studies show that if prompt explo- 
sions occur, one objection to prompt r-process nucle- 
osynthesis is avoided. There is no over-production of 
r-process elements as long as an accurate accounting 
of the neutronization of the prompt ejecta is made. 
It also provides a very good reproduction of the Solar 
r-process abundance distribution for elements with 
A > 100, suggesting that this might be the source 
of the heaviest r-process nuclei while perhaps delayed 
neutrino-driven winds or neutron-star mergers may 
be the source of lighter r-process nuclei. Obviously, 
a key component of this paradigm is to understand 
just which progenitors are capable of producing a 
prompt explosion and r-process nucleosynthesis and 
also whether the late time neutrino heating will affect 
the r-process yields. In future work we will address 
both of these points, however, regarding the delayed 
neutrino emission we expect that it will not affect the 
results reported here as by the time the delayed neu- 
trinos arrive, the prompt ejecta is already far from 
the proto-neutron star. 
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Fig. 2. — Evolution of: (a) baryonic mass density; (b) temperature; (c) entropy; and (d) electron fraction as a 
function of time for ejected material in the prompt explosion. Note that the time interval is expanded for (d) to 
better show the rapid change in Yg. 
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3. — Electron fraction as a function of baryon mass coordinate for ejected material in the prompt explosion. 
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Fig. 5. — Final integrated isotopic abundances for ejected material (solid line) compared with the Solar r-process 
abundances (filled circles; Kappeler, Beer, fc Wisshak 1989| ). 
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